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We present data for the scaling behavior of lattice QCD with two flavors of light Wilson fermions. This is done 
by matching 7r n and p masses at the two lattice sizes of 16 x 32 and 8 x 16. We find that at /3 = 6/g = 5.3 
the matching does not occur over a range extending down to /3 = 3.5. For f3 = 5.5 matching may be achieved at 
f3 = 4.8 — 4.9, leading to a A/3 = 0.6 — 0.7 which is higher than the perturbative 2-loop value of 0.45. In both 
cases we conclude that the simulations are very far from the perturbative scaling region. 



1. Introduction 

If we measure some physical observable fi, such 
as a hadron mass, in a Monte Carlo computa- 
tion of lattice QCD the result will not depend on 
the details of the lattice discretization as long 
as the lattice spacing a is small enough. This 
is just the assumption that there is a universal 
continuum limit of lattice QCD. The value of 
the observable obtained by numerical calculation 
fi = a s f(g,K) depends only on the parameters 
appearing in the lattice action, the gauge cou- 
pling 6 = 6/g 2 and the hopping parameter re. 
The dimensionless quantity / has been related to 
the dimensionful observable fi by multiplying by 
the appropriate power of the lattice spacing a. As 
we approach the continuum limit a — > the fact 
that fi is independent of a is expressed by the 
renormalization group equation 
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As we have two independent relevant parame- 
ters we require two independent physical observ- 
ables to use as renormalization conditions to fix 
these parameters. 

In a perturbative expansion of /3i about g = 
only the first two terms are independent of the 
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details of the lattice regularization scheme. 

The aim of this work is to determine how far the 
values of g being used are from those at which /3i 
is perturbative for the realistic case of QCD with 
two flavors of light Wilson fermions. 

If we select some set of physical observables fi^ 
we have 
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where 8i is the dimension of fi^ and fi is the corre- 
sponding dimensionless function measured on the 
lattice. In the continuum limit all these observ- 
ables must become independent of the lattice, so 
for two different lattice spacings oti and a-i 
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In practice we work with two lattices, Ai with 
L\ sites in each direction and with lattice spacing 
oi, and A2 with L 2 sites in each direction and 
lattice spacing a 2 . By fiat we assert that Ai and 
A2 represent the same physical volume, thus 



L\a\ = L 2 a 2 . 



(5) 



Using this relation we recast Equation (4) in 
terms of dimensionless quantities only 
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We define the quantity 
A/3: 6 6 
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to be change in the coupling needed to compen- 
sate for this change in the cutoff. 

2. The /3-function 

The /3-function at fixed re 1 is defined by f3i(g) = 
—adg/da, and the two universal terms from per- 
turbation theory are f3i(g) = —bog 3 — big 5 + • • •. 
The solution of these equations, expressed in 
terms of /3 = 6/[g(a 2 )] 2 , is 
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after evaluating the integral we find that 
A/3 Peii 
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which is the perturbative prediction for the 
change in (3 needed to effect a change in scale 
by a factor aila\. 

3. Previous Results 

The (3 function for two flavors of Wilson 
fermions was computed by some of us earlier [1]. 
The physical observables used for that purpose 
were a 7r-like state and several Creutz ratios of 
Wilson loops. 

The main conclusion then was that perturba- 
tive scaling is expected to occur only at (3 > 6.0. 
In those simulations the 7r-like object had a heavy 
mass and it was thought at the time that this may 
have hidden the effects of the dynamical fermions 
in the system. Simulations with lighter 7r masses 
are thus necessary to clarify the issue. 

4. Present Simulations 

In this work we concentrated on data available 
on lattices of size 16 3 x 32 with two flavors of Wil- 
son fermions with parameters and results shown 
in Table 1. 

The values of (3 were chosen for hadron spec- 
trum and other calculations largely based on the 
assumption that the presence of two fermion fla- 
vors renormalizes /3 que nched downwards by about 



We expect the /3-function to become independent of k in 
the continuum and chiral limits. 



Table 1 

Parameters and Masses on 16 3 x 32. 



5.3 0.1670 0.4540(20) 
5.3 0.1675 0.3120(40) 
5.5 0.1604 0.2615(31) 



am p m^/mp 

0.6350(20) 0.7150(54) 

0.5230(40) 0.597(15) 

0.4278(70) 0.611(20) 
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Figure 1. Crosses: 
(16 3 x 32). 



(3 = 3.5 (8 3 x 16); □: f3 = 5.3 



0.5 and thus we may not be far from the pertur- 
bative scaling regime at these parameter values. 
For our test of scaling we used two hadron 
masses (m w and m p ) as observables for match- 
ing with similar data generated on lattices of size 
8 3 x 16 at values of /3 ranging from 5.1 down to 
3.5 and many re values. 

5. Results 

Our first result is the unmatchability of the 
data at (3 = 5.3. We find here that A/3 must 
be larger than 1.8. Figure 1 shows the two points 
obtained at (3 = 5.3 on the larger lattice (with 
m w appropriately scaled by a factor of 2) and the 
corresponding data for m w and m^/mp on the 
smaller volume for (3 = 3.5. 

Our second result is that preliminary analysis 
of our data shows that at (3 = 5.5, matching is 
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Figure 2. *: (3 = 5.5 (16 3 x 32); o: (3 = 4.8, 
o: /3 = 4.9 (8 3 x 16). 
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possible (Figure 2) with points at (3 between 4.8 
and 4.9 on the smaller lattice. We find for this 
case therefore that A/3 = 0.6 to 0.7. 

A third result is an approximate parameter re- 
duction in the theory at strong coupling: (3 = 3.5 
to (3 = 4.7. Figure 3 includes all our data in this 
range of (3 on the smaller lattice and shows that 
for 3.5 < (3 < 4.9, the data is essentially a func- 
tion of a single parameter. A change in (3 may be 
compensated for by a change in re. 

6. Conclusions 

Combining the present results with those ob- 
tained earlier on smaller lattices one is led to the 
conclusion that current dynamical simulations are 
being done at strong coupling. This is clear from 
the /3-function we have measured, which is shown 
in Figure 4. 

One therefore may need to reassess estimates 
for parameters and hence the cost of future com- 
putations with dynamical Wilson fermions ac- 
cordingly. 

Acknowledgements 

This research was supported by by the U.S. De- 
partment of Energy through Contract Nos. DE- 
FG05-92ER40742 and DE-FC05-85ER250000. 



Figure 3. Cummulative data for 3.5 < (3 < 4.7 
(8 3 x 16). 




Figure 4. QCD /3-function for m^/mp ~ 0.6 (o) 
and from Ref. 1 (—). 



REFERENCES 

1. K. Bitar, A. D. Kennedy and P. Rossi, 
Phys. Rev. Lett, 63:25, (1989), 2713. 



